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We study theoretically 2D HgTe/CdTe quantum well topological insulator (TI) illuminated by 
circularly polarized light with frequencies higher than the difference between the equilibrium Fermi 
level and the bottom of the conduction band (THz range). We show that electron-hole asymmetry 
results in spin-dependent electric dipole transitions between edge and bulk states, and we predict an 
occurrence of a circular photocurrent. If the edge state is tunnel-coupled to a conductor, then the 
photocurrent can be detected by measuring an electromotive force (EMF) in the conductor, which 
is proportional to the photocurrent. 


I. INTRODUCTION 

Topological insulators (TIs) became a focus of atten¬ 
tion of many condensed matter physicists in recent years, 
not least due to their possible applications in spintron- 
ics and quantum computing. These are materials with 
time-reversal symmetry and non-trivial topological or¬ 
der, which have an insulating bulk but conducting topo¬ 
logically protected edge/surface statei^H^. Spin-orbit in¬ 
teraction plays a significant role in these materials, and 
particularly manifests itself in spin-momentum locking of 
charge carriers in edge/surface states. 

Optical excitation is an efficient tool for generating cur¬ 
rents in materials. This process has been studied in 3D 
TI^^Hm In Ref. m effects in study are due to electric opti¬ 
cal transitions between surface states, which are possible 
in the presence of magnetic field. Besides, in the presence 
of strong magnetic fields, the electric di^le transitions 
between Landau levels are also possiblei^. In contrast 
to these papers, here we study 2D TIs and show that 
the optical generation of the current is possible without 
magnetic field as well. 

HgTe/CdTe quantum well structurei^^ are one of the 
most well-known 2D TIs. These quantum wells exhibit an 
inverted band structure if their width exceeds a certain 
critical value. The inverted band structure and strong 
spin-orbit interaction give rise to unusual optoelectronic 
phenomena, e.g. a nonlinear magneto-gyrotropic photo- 
galvanic effect (PGE)P^ (PGE). A circular PGE was also 
experimentally observecpSl when the sample was illumi¬ 
nated by mid-infrared or terahertz laser radiation. These 
photocurrents were induced due to direct transitions be¬ 
tween different size-quantized subbands or due to indirect 
(Drude-like) transitions within the lowest size-quantized 
subband. In both cases the optical transitions responsi¬ 
ble for the PGE involve only bulk states. However, in 
the case of a finite sample size of 2D TI there exist topo¬ 
logically protected helical edge states which form two 
branches with opposite spin^i^. In the recent papei^ 
it was predicted that edge states affect bulk magneto¬ 
conductivity. The paper focuses on the bulk properties 
in a strong magnetic field and does not discuss opto¬ 


electronic properties of the edge states. However, it is 
of interest whether the photocurrent can be induced at 
the edge states at zero magnetic field. A PGE due to 
transitions between edge states of the opposite chiralities 
has been predicted in Ref. [5D1 but electric dipole transi¬ 
tions between them are forbidden by selection rules, and 
only magnetic dipole transitions are possible in this case. 
Thus, the direct transitions between the edge states are 
weak. Unlike Ref. HOI in this paper we study electric 
dipole transitions between the edge and bulk states in 
HgTe/GdTe quantum well 2D TI in zero magnetic field, 
which lead to the edge currents. To our knowledge, this 
mechanism of the PGE in TI has not been studied yet. 

Starting from the Bernevig-Hughes-Zhang (BHZ) 
model of HgTe/GdTe 2D TI, we find a relation between 
matrix elements of the edge-bulk transitions. In the gen¬ 
eral case the electron-hole symmetry is broken, and the 
probability of transition depends on the spin and, hence, 
on the chirality. Thus, the transitions will lead to a dif¬ 
ferent population of spin-up and spin-down states and 
to occurrence of a photoinduced electric current. In or¬ 
der to study this effect, we derive a kinetic equation and 
then solve it in the quasi-equilibrium approximation. We 
also propose a way to detect the photoinduced current by 
coupling 2D TI with a conductor and measuring an elec¬ 
tromotive force (EME) induced in the conductor. 

The paper is organized as follows. In Section|TT]we con¬ 
sider optical transitions between edge and bulk states and 
derive the photoinduced electric current. In section [llT| we 
calculate the EMF that appears in the conductor tunnel- 
coupled to the edge state. 

Below we set h = 1, c = 1. 

II. OPTICAL TRANSITIONS BETWEEN EDGE 
AND BULK STATES 

We consider a HgTe/CdTe quantum well TI with a 
conducting helical edge state illuminated by circularly 
polarized light with a frequency ujq slightly exceeding the 
the absorption threshold, so that optical transitions may 
occur between the edge state and the bulk conduction 
band (Fig. [^. We assume that the TI is located at a: > 
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FIG. 1. Schematic picture of optical transitions between two 
edge state branches and bulk states. Circularly polarized light 
of frequency ujo exceeding the absorption threshold by Au 
induces transitions between edge and bulk conduction band 
states. The different thickness of arrows depicts that probabil¬ 
ity of electric dipole transitions may depend on chirality/spin 
of the edge electron. 


0; j/-axis is taken along the edge of TI, and z-axis is 
perpendicular to the 2D TI. 

Both in HgTe and CdTe the relevant bands are s-type 
band (Tg) and p-type band split by spin-orbit interaction 
into a J = 3/2 (Tg) and a J = 1/2 (Tr) bands. The lat¬ 
ter is usually neglec ted as it has negligible effects on the 
band struct ure^^I^. CdTe has a band-ordering similar to 
GaAs with Tg conduction band and Tg valence band. In 
HgTe the usual band-ordering is inverted. The quantum 
well subbands derived from the heavy-hole Tg are usu¬ 
ally denoted by i?„, and the subbands derived from the 
electron Tg are denoted by 

We describe HgTe/CdTe quantum well TI by the four- 
band BHZ modeP^. 

In the four-component basis consisting of \Ei f), 
\Ei i), \H,i) with mj = 1/2,3/2,-1/2,-3/2 
correspondingly, the Hamiltonian reads: 

Hbhz = —T’k^ + 

Ak+ 0 ° \ 

Ak_ 0 0 

0 0 M-Bk"^ -Ak_ 

Vo 0 -Ak+ -{M-B\A)) 

( 1 ) 

where k± = k^ A iky. Here A, B, V, M are material 
parameters, which depend on quantum well geometry: 
A > Q, B < Q] parameter M. is negative if the quantum 
well is in a TI state, and 2|Ad| is a value of a band gap in 
TI. Parameter I? 0 if electron-hole symmetry is broken, 
and D = 0 otherwise. 


Since the Hamiltonian Q has been obtained in kp- 
approximation, and the non-diagonal part of Q linear by 
quasi-momentum k corresponds to kp-term, the Hamil¬ 
tonian of a light-matter interaction in electric dipole ap¬ 
proximation reads 

He-A = Hi {Aj: — iAy) -|- hI [Ax iAy) (2) 

Hi = |e| \Hi t) {El t| - |e| \Ei i) {Hi /| (3) 

where e is the electron charge, A is a vector-potential of 
electromagnetic field. In the case of right-hand polarized 
light (as defined from the point of view of the source) 
propagating along z-axis the vector-potential can be rep¬ 
resented as 

= -cosa;oc, Ay = -sincjQC, A^ = (J, 

UrOJo TlriOo 

( 4 ) 

where W is the intensity of light, Ur is the refractive 
index. 

Under the illumination by the right-hand circularly po¬ 
larized light the selection rules allow only those electric 
dipole transitions from the edge to bulk, which increase 
angular momentum by Tl. In the Hamiltonian ([^ these 
transitions are described by the first term. The conjugate 
term describes the reverse optical transitions. 

In the absence of the boundary, the eigen states of 
2D TI are bulk states separated by a gap, and the bulk 
states at the bottom of the conduction band are formed 
by Hi states with a zero momentum and a well-defined 
projection of an angular momentum mj = ±3/2. 

In a finite-size sample, there appear edge states which 
are superpositions of Ei and Hi Bloch wavefunctions. 
These edge states can be found by solving Schroedinger 
equation with zero boundary conditions for the Hamilto¬ 
nian Q in the coordinate representation in ^-direction 
and the momentum representation in y-direction (see 
Ref.HH): 

ij’edgeA OC (^iSl/lB-V] \Ei,s) + l/lB + V] |i?l, s)^ X 

(e-U,-" - e-^-+") , (5) 

where are inverse decay lengths for the localized 
edge states. 

In the presence of the boundary, conduction band 
wavefunctions are distorted near the boundary where 
they overlap with the edge states. Instead of explicitly 
calculating the conduction band wavefunctions by solving 
Schroedinger equation with zero boundary conditions, we 
will use the time-reversal symmetry and orthogonality 
conditions for the eigen states of the Hamiltonian. 

The BHZ Hamiltonian and zero boundary condi¬ 
tion are invariant under the time-reversal symmetry 0. 
Therefore, if the spin-up eigen state with energy e is of 
the form ky,x) = fe,ky{x) \E t)+ge,ky{x) \H f), 

then the spin-down eigen state can be obtained as 

Q'^buik.t = -fl-ky{x) \Ei j.) - gl_ky{x) \Hi j,) 
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We denote the overlap integrals between the edge and 
bulk states as: 

+ 00 

F{e, ky) = j - e~^+^) {x)dx, 

0 

+ 00 

G{^,ky)= J g^^ky{x)dx. 

0 


Mutual orthogonality of edge and bulk states yields the 
relation between the integrals 


F_ ^ yig+p| 

G 


( 6 ) 


Thus, both / and g are non-zero, and not only the edge 
states are superpositions of E and H Bloch wavefunc- 
tions with different well-defined projections of total an¬ 
gular momentum, but the bulk conduction states are also 
their superpositions even at the bottom of the band. Se¬ 
lection rules allow transitions from |i? t) to t)i and 
from \H j,) to \E 1). Thus, the transitions from the both 
edge state branches to the conduction band satisfy the 
selection rules. Calculation of the matrix elements Wg of 
the first term of the Hamiltonian ([^ corresponding to 
these transitions yields the main relation 


K| ,/\B^\\G\ _B-V 
ktl ^\B + V\\F\ B + V 


of order of 1 meV and V jB is order of 0.5. These values 
of parameters correspond to characteristic frequencies ojq 
of order of 100 GHz and ratio of matrix elements 
of order of 10, i.e. the probabilities of transitions is also 
strongly spin-dependent. 

Note, that we used zero boundary conditions (BCs) 
for 2D TI. The result Q does not qualitatively depend 
on the BCs for the wavefunctions provided they are in¬ 
variant under time-reversal symmetry and yield helical 
edge states. Different BCs are discussed in Ref. EHETl 
Although the choice of the boundary conditions does not 
affect the topological nature and the existence of the edge 
states, the spectra of bulk and edge states, and the eigen¬ 
states themselves depend on the BCs. Particularly, the 
matrix elements of the transitions may depend on the 
BCs. However, in our approach only the Eq. ([^ depends 
on the BCs. In general case the amplitudes of Ei and JIi 
Bloch functions will be different, but if the electron-hole 
symmetry is broken these amplitudes will remain still un¬ 
equal. Further, in order to obtain our main result 0 we 
exploit their inequality, time-reversal symmetry and mu¬ 
tual orthogonality of eigenstates. Thus, we believe that 
the result does not qualitatively depend on the BCs pro¬ 
vided they are invariant under time-reversal symmetry 
and yield helical edge states. 

If the light is incident in an arbitrary direction n®’'^ = 
(cos(?!)sin0,sin(/)sin0,cos0), then the matrix elements 
wfcan be obtained by replacing vector potential in the 
light-matter interaction Hamiltonian ([^ with its projec¬ 
tion on the TI plane(for details see Appendix [C|): 


Note, that if electron-hole symmetry is present, the prob¬ 
abilities of optical transitions from the both edge state 
branches are equal. Electron-hole symmetry implies that 
the edge states are superpositions of \E,s) and \H,s) 
with equal (up to a phase factor) amplitudes. The same 
is true for the bulk states (it is shown explicitly in the 
Appendix 0. Hence, the transition probability for the 
both spin-up and spin-down branches will be the same. 
However, in real samples the electron-hole symmetry is 
broken i.e. the electron and hole components of the bulk 
eigen states as well as the edge eigen states are not equal 
anymore, and the probability of an edge-bulk transition 
will depend on spin. 

Values of matrix elements for the case of strong 
electron-hole asymmetry are derived in Appendix [Bj 

We estimate the ratio of probabilities for typical val¬ 
ues of parameter^ A = 365 meV • nm, B = —686 meV • 
nm^,I? = — 512 meV • nm^, Ai = —10 meV correspond¬ 
ing to the quantum well width dc = 7 nm. In this 
case « 47.4. Thus, the transitions in case 

of HgTe/CdTe 2D TI are strongly spin-dependent. 

Our consideration can be applicable not only in case of 
HgTe/CdTe 2D TI but also in case of other 2D TI which 
can effectively described by BHZ model. One of the in¬ 
teresting examples is a recently predicted all electron TI 
in InAs double welP^ which allow easily tune BHZ pa¬ 
rameters. In case of this material the band gap 2|A4| is 


,e,4 >_Wse *'^(I-fcos6»)-Pw*^e*'^(I-cos6<) 

'2 ' ' 

Below we will use an expression for the value of the 
squared matrix element averaged over the direction: 

= + ( 9 ) 

Kinetic equations for distribution functions of electrons 
in the edge state n{e) and in the conduction band N{e) 
can be written as 

dnsje) ^ ns{e) - Ng{s + ujq) ^ 
dt Tind,s{e,£ + Xlo) 

"■("+") [ 1 -’;■(")! A.+ 

Tgp{£ + UJ,£) 

n-sje) - ns{£) 

Te 

where Wr//); is the rate of transitions induced by illumi¬ 
nation, Tg“^(e',e) is the rate of spontaneous transitions 
between the conduction bulk state with energy e' and the 
edge state with energy e, i.e. recombination rate; is the 
spin relaxation time for the edge electrons. Since there 
is still discussion in the literatur^^MlIl which mechanism 
agrees better with experimental data, in this paper we 
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introduce this time assuming that in any realistic system 
it is finite. 

Here we write these kinetic equations ( |10[ ) phenomeno¬ 
logically, and more rigorous derivation based on the 
Keldysh technique is given in the Appendix The rate 
of the transitions induced by illumination can be related 
to the matrix elements using the Fermi golden rule 

where vc,ky is the density of states in the conduction 
band with a fixed ky-. 


VC,k, 


~ ^ ^ ^k^,ky) — 




\/2-K^e — \M\ — fc^/(2m) 


where m is the effective mass of conduction band elec¬ 
trons and Lx is the length of 2D TI in the x direction. 
Note that summation over kx in the definition of i'c,k 
instead of summation over both kx and ky reflect the 
fact, that the transitions are vertical, i.e. ky projection 
of momentum conserves. 

In order to deduce a relation between induced and 
spontaneous transition rates one can use a detailed bal¬ 
ancing condition similar to that for Einstein coefficients 
for discrete level^^. The factors Tgp^ind in the kinetic 
equation do not depend on the illumination and environ¬ 
ment, since they are intrinsic properties of the 2D TI. 
Therefore, the kinetic equation (10) should remain valid 


if we put the system in thermal equilibrium with black 
body radiation. In this case the distribution functions 
of the edge and bulk electrons are the equilibrium Fermi 
function with the same Fermi level, and photons have 
the Bose distribution. The detailed balancing between 
the states with energy e and e' = e -I- w for an arbitrary 
ui yields 


Finally, we obtain the exp ress ion for the spontaneous il¬ 


lumination rate from (11) 


Tsp^ = 3 + k-sH nruji>c,ky{e + uj) (12) 


We solve kinetic equation (10) in the quasi-equilibrium 


approximation, i.e. assuming that the distribution func¬ 
tion of edge electrons with spin s is a Fermi distribution 
with a quasi Fermi levels sp + fJ-s-, and, similarly, the dis¬ 
tribution function of conduction bulk electrons with spin 
s is a Fermi distribution with a quasi Fermi level |A1| -l-Cs 
(|A1| is the bottom of the conduction band). The quasi¬ 
equilibrium approximation can be justified if the energy 
relaxation times in the edge and bulk states are much 
shorter than the life-time of the excess photogenerated 
electrons. Since the results will depend on whether the 
initial Fermi level is above or below the Dirac point we 
consider both these cases. 


A. Fermi level above Dirac point: absorption 
without photocurrent 


The electrons in the conduction band in the quasi¬ 
equilibrium approximation lie in the bottom of the con¬ 
duction band and they can recombinate only with empty 
states in the vicinity of the Dirac point. If the Fermi 
level is above the Dirac point (see Fig. [^) then all the 
states near Dirac point are occupied 1 — ns(e) = 0, and 
spontaneous transition term in (llOl) vanishes. 


After integrating the kinetic equation (10) over ener¬ 


gies we obtain a relation between quasi-Fermi levels 




2,.a 


ntuj. 


r^O 


[no(£) - no(£ -f w)] {Tind,s)g^^ 

[no(£) - no(£ -f W)] {nnd,-s)g^^ ^ 

no{s + uj) [l-no(£)] 

'^sp 

where is the induced transition rate averaged 

over the direction of an incident equilibrium photon, 
no = ^ (l ~ tanh is the equilibrium Fermi distribu¬ 
tion, and = ^^^Nph{u}) is the spectral density 

of equilibrium right(left) polarized illumination with the 
photon distribution function Nph{u)) = | (coth ^ — l). 
The value of the induced transition rate (r^”^ averaged 
over direction can be calculated using 

_ 1 \ _ ind,s ' ind, — s 

’^ind) g,rt> 


where Aw is the difference between the light frequency 
and the absorption threshold Aw = ep + ooo — \M.\ — 
ky/{2m) (see Fig. flj) The only solution is d-s = ds = 
—Aw. Almost affi the electrons with energies from 
Ep — Aw to £i:’ are moved to the conduction band by 
illumination(/r^ « —Aw). They equilibrate in the con¬ 
duction band and remain there, since the edge states 
near the Dirac point are occupied. The electric current 
j = Go(/r^ — dl) = 0) where Gq = jh is the conduc¬ 
tance quantum. 


B. Fermi level below Dirac point: non-zero 
photocnrrent 

The situation differs if the Fermi level is below the 
Dirac point. In this case the edge states in the vicinity 
of the Dirac point are not occupied and the spontaneous 
transitions from the conduction band to edge states are 
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FIG. 2. (a) Schematic picture of quasi-Fermi levels if Fermi level is above the Dirac point. All the spin up and spin down 
edge electrons with energies from ef — Ao; to ef are moved to the conduction band. States near the Dirac point are occupied, 
and, hence, spontaneous transitions from bulk to edge states are not allowed. The electric current of the edge electrons 
j = Go(Ai-t- — = 0; (b) Schematic picture of quasi-Fermi levels if Fermi level is below the Dirac point and for low intensities 

VF < Wj. — Wt (region I on the Fig. |^. All the spin-up electrons with energies from ef — Atj to ef are excited by the light, and 
spin-down electrons in the conduction band appear mainly due to spin relaxation. Recombination of bulk spin-down electrons 
shifts the quasi-Fermi level of spin-down electrons above the initial Fermi level, (c) Schematic picture of quasi-Fermi levels for 
high intensities W > VV’^/2 (region II on the Fig|^. Almost all spin-up and spin-down electrons with energies from e — Aw 
to E are excited by the light(/i.|. « rs —Aw). The spin imbalance at the edge states decays with intensity of the light as 
Sfi oc W~^ 


allowed. Integration of (10) yields 


yv^/JI^+~Aw = fc + To— —(13) 

V AcJ \ '^e / 


m = 


Sttfc 


Wo i\M\-eD) 


UTi 


Aw 


ksp + 3 


8L^nrmvTi{\M\-eD)(\w\'^s + \M-s) 

■ 


where Aw = ep + wq — |AI| (see Fig. [^. Here I'e and 
VC are the densities of the edge states and the conduc¬ 
tion bulk states correspondingly. Note that in contrast 
to vc^ky defined above, vq ^C,ky is the 2D density 

ky 

of states in which summation over both components of 
momentum is performed. Another relation results from 
the conservation law for the number of particles 


^e(Mt + + Cf) — 0- (14) 


In the stationary regime neither spin nor charge accumu¬ 
lates, and the same number of transitions per unit time 
occurs from spin-up states to spin-down states and vice 
versa. Therefore we can equate spin-relaxation rates of 
the edge and bulk electrons and obtain another quasi¬ 
equilibrium condition: 


Ve 


Ms M—s 


Te 


VC 


C-s - Cs 
TC 


(15) 


where tc ^ is the spin relaxation time for the con¬ 
duction band electrons. 



FIG. 3. Dependence of photocurrent (solid line) and currents 
of spin-up/spin-down electrons (dashed lines) on light inten¬ 
sity if the initial Fermi level is below the Dirac point. The 
currents of each branch are proportional to the corresponding 
quasi-Fermi level shift. For the region I (W < 1 ^ 4 ./2) edge 
branches contribute to net current with the same sign, while 
in the region II (W > VV’^/2) the current of spin-down elec¬ 
trons changes its sign. Currents are measured in the units of 
Go Aw. 


An important limiting case is when one can neglect 
spin relaxation of the edge electrons assuming that Tg ^ 
Tc a nd bot h times are great enough. In this case equa¬ 
tion 


14 


15 


yield = C-s = equa- 
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tion (13) can be solved analytically. 


^^o — 






2W2 (w^ 




(Wf 


wf) 


Au! 

(16) 

(17) 


. W|-W| 2 

An electric current j = GqS^j, arises in the edge state, 
where Gq = jh is the conductance quantum. If the in¬ 
tensity of light is small W <C (Fig[^), mainly the 

spin-up electrons are excited by the light. The spin re¬ 
laxation time for the bulk states is much shorter than the 
life-time of photogenerated electrons, hence, spin of the 
photogenerated electrons relaxes, and there appear addi¬ 
tional spin-down electrons in the conduction band, whose 
recombination rate turns out to be greater than the exci¬ 
tation rate of spin-down edge electrons. Thus, the quasi 
Fermi level for spin-up electrons is below eF^ and the 
quasi Fermi level for spin-down electrons is above it. The 
currents of electrons with the opposite spin contribute 
to the total electric current with the same sign (see 

Fig Sr, Figj^), and the total current j « 2Gq Auj. 

I: the intensity of light W » (Fig.[^), all the 

edge electrons with energies from \M\—ujq to ep are ex¬ 
cited to the conduction band, and the quasi Fermi levels 
of the edge electrons saturate at the value = —Aw. 
The currents js = Gofig are almost equal, but contribute 
to the total current with opposite signs. Thus, the to¬ 
tal current decreases with the increase of the intensity as 
j oc }V~^. For an HgTe/CdTe quantum well with width 
dc = 7 nm, and the sample size of order L ^ 1 /rm, 
the quasi Fermi level for spin-up electrons saturates at 
W ~ ~ 10“® W/cm^. The quasi Fermi level for 

spin-down electrons saturates at ~ 10“® W/cm^ 

However, in case of extremely weak intensities spin re¬ 
laxation rate turns out to be comparable with transitions 


rates and hence cannot be neglected. We can solve (131- 


(15) analytically assuming W < ‘C and obtain 


j = Go Auj 


tI 

4W2 ^2 


1 - 


/ Wt 


The full curve for the dependence of photocurrent on 
intensity is sketched on the Fig. 


III. TUNNEL CONTACT TO AN EXTERNAL 
CIRCUIT 

In the previous section we showed that circularly polar¬ 
ized illumination induces an electric current at the edge of 
HgTe/CdTe quantum well TI. In order to observe the ef¬ 
fect one should connect the sample to an external circuit. 
We consider a system of HgTe/CdTe tunnel-coupled by a 
contact of length A to a 2D metal conductor (see Fig. |^. 



FIG. 4. A proposal to detect photocurrent in the edge state; 
2D conductor of length Ly is coupled to the edge state via 
the tunnel contact of length A. The EMF appears between 
the opposite ends of the 2D conductor. 


The Hamiltonian of the system reads 

H = Htj + H^d + Htuu (18) 

In order to describe the edge state we use an effective 
edge state Hamiltonian corresponding to the linear edge 
state spectrum 

Hti = X! / '^l^y)i-™Tisdy + eD)ip8{y)dy, (19) 


where ip is the effective field operator for an electron in 
the edge states, is the energy of Dirac point mea¬ 
sured from the middle of the band gap, vpi is the ve¬ 
locity of edge electrons. We assume that illumination 
results in quasi-equilibrium occupation numbers of elec¬ 
trons corresponding to the quasi Fermi levels, and the 
tunnel coupling is weak that it does not affect the opti¬ 
cal transitions, so the results of the previous section are 
applicable. The Hamiltonian of the 2D conductor reads 

H2D = H2D,o + ^ y (20) 


where H 2 d,o is the Hamiltonian of free 2D electrons, df is 
the field operator in the 2D conductor. Here we take into 
account a random delta-correlated potential of impurities 
Vimp characterized by a mean scattering time Timp- 
The tunneling Hamiltonian reads: 

Htun = j dy ipl{y)T'^six = 0 ,y) + H.c. (21) 


where T is a matrix element of tunneling, and we assume 
that the tunneling is momentum-conserving. 

We start from the Hamiltonian (18)-(21), and then 
derive equations for Keldysh Green functions 
taking into account the tunneling Hamiltonian and the 
impurity potential as perturbations. The self energy of 
2D electrons in the conductor resulting form the tunnel¬ 
ing reads 


S^^,f^’^(i‘i,r2) = |Tp(5(a;i)5(a;2)Gyj^^’'^(yi - j/2) 
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The Green functions of the electrons in the edge state 
electrons can be obtained by solving the corresponding 
Dyson equations: 

GT\^]{£,Py) = 

Grpsi^^Py) = 

(23) 


' ( 22 ) 
_ (^-^1 _ 

-27rt(5r(e - £p^J tanh ^- — 


where T is the inverse life-time of electrons in the edge 
state, which consists of contributions from tunneling 
and optical transitions. The contribution from the op¬ 
tical transitions is determ ined b y self energy operator 
iT = J2kyi'^^ ~ ^^) (D15), (D17|). Besides, some 

mechanisms of spin-relaxations, e.g. coupling to multi¬ 
ple puddle^^ also contribute to T, so here we introduce 
it phenomenologically. Note, that the exact value of T 
does not affect the final results. The Lorenz-type factor 
i5r {x) = ^ y. 2 ^r 2 describes broadening of electronic states 
due to their finite life-time. 


The derivation of kinetic equation for quasi-classical 
distribution function /s(r,p) in the 2D conductor is 
straightforwarcP^ and gives 


[dt + (v, V)] fs+ 

2'K\T\'^6{x)e - |y|^ ,5r [fs - Us] = 


Here Us is a distribution function for the edge state elec¬ 
trons the TI, and the Heaviside step-function 6 restricts 
the length of the contact to A. 

If the TI is illuminated by circular-polarized light, the 
quasi Fermi level of spin-up electrons is lower than the 
quasi Fermi level of spin-down electrons. In a stationary 
regime the current through the tunnel contact should be 
equal to zero, therefore, the Fermi level in the 2D con¬ 
ductor is exactly in the middle between the quasi Fermi 
levels in the edge states of TI. However, the zero tun¬ 
nel current consists of a spin-up current from the TI to 
the conductor and a spin-down current from the conduc¬ 
tor to TI. Electrons of the opposite spins in the TI are 
of the opposite chiralities and the tunneling is assumed 
to conserve momentum,so it should result in an electron 
drift and appearance of the counterbalancing EMF in 
the conductor. We consider a stationary case in which 
tunneling leads to charge redistribution, and, hence, the 
extra charge induces an electric field.The electric field E 
is related with the extra charge —\e\5p by the Poisson 
equation: 


div E = —47r|e|i5p, (25) 

where Sp is the deviation of electron density from its 
equilibrium value, and the full electron density can be 


expressed in terms of the distribution function as 



The kinetic equation that takes into account the elec¬ 
tric field in the conductor reads 

[5t + (v,V)]/« + n, + ^^^^-|e|vE^ =0, (27) 

Hs = 0 - lyl^ 2Tr\T\'^S{x)6r (sp^^, - [fs - 

(28) 

Here Hg denotes the term responsible for the tunneling. 
Since the electric field E arises due to the deviation of the 
distribution function / from its equilibrium unperturbed 
value /o, we can replace fs in the field term with /o, 
neglecting second-order corrections. 

We consider a model where the 2D conductor is a nar¬ 
row strip spanning from x = —Lxl2 and x = and 

is infinite in y-direction, and the contact is in the middle 
of the strip at x = 0. For simplicity we assume that the 
width Lj, of the conductor is smaller than the scattering 
length vpT. In a more general case the conductor can be 
qualitatively considered as the same narrow strip with 
a bypass resistance. Under these assumptions one can 
take into account only the electric field along the contact 

Ey » Ey,. 

We solve the kinetic equation together with the Pois¬ 
son equation assuming quasi-neutrality and treating the 
angle-averaged distribution function fg in the impurity 
term self-consistently (for details see Appendix . 

Finally, we obtain the EMF as the integral over the 
conductor length of the electric field averaged in the 

ix/2 

transversal direction {Ey)^ = L~^ f Ey dx: 

-ix/2 


Ly j‘1 

\e\E= J \e\{Ey)^dy = 

-Ly/2 

irp A Sp 

where pp, vp are the Fermi momentum and velocity in 
the conductor, and 6p is the difference between quasi 
Fermi levels given by 0 

IV. CONCLUSIONS 

In summary, we considered electric dipole optical tran¬ 
sitions from helical edge states of HgTe/CdTe TI illumi¬ 
nated by circularly polarized light to the bulk conduction 
band. If electron-hole symmetry is broken which typ¬ 
ically is the case, these optical transitions are strongly 
spin-dependent. 











This gives rise to a circular electric current in the edge 
state if the Fermi level is below the Dirac point. The 
value of the photocurrent reaches maximum and then 
decreases with the growth of the light intensity. It is 
worth noting that although the overlap between wave- 
functions of the edge and bulk states determines the time 
required for stationary regime to settle in, the magni¬ 
tude of the current in this regime does not depend on 
the overlap, and, hence, we anticipate that the effect can 
be observed even in samples where this overlap is small. 
We showed that the photocurrent can be detected elec¬ 
trically by measuring EMF in the conductor coupled to 
the edge state of the TI. 
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Appendix A: Optical transitions in the case of 
electron-hole symmetry 


of incident and reflected waves, while the second terms 
correspond to the parts of wavefunctions localized near 
the boundary. In the presence of the boundary the angu¬ 
lar momentum is ill-defined. Although the wavefunctions 
of the conduction band behave like \Hi,mj = 3/2) and 
\Hi,mj = —3/2) with a well-defined angular momentum 
far away from the boundary at x ^ k~^, the overlap 
integral between edge and bulk states is dominated by 
small distances from the boundary x ~ AZ^, where these 
wavefunctions are superpositions of \Ei), \Hi) with equal 
(up to a phase factor) amplitudes: 




i 

0 

Voy 

0 

0 

1 

—V 


(1 - 




(A4) 


(AS) 


Thus, the matrix elements are equal up to a phase factor 

ktP = kfl" ~ L^- 


In this section the matrix elements of the optical transi¬ 
tions are explicitly calculated for the case of electron-hole 
symmetry V = 0. In this case, the Hamiltonian Q with 
zero boundary conditions yields the edge state wavefunc- 
tion 


'^edge,s{^: ky') — edge^s (e C ^ , 


edge,'[ ^ 



"^edge.X OC 



(Al) 


where A+ « |^, A_ « <C A+. Here we assumed that 

^ BA4, which is the case for the typical parameters 
of TI. Solving Schroedinger equation with zero boundary 
condition for e — |M| ^ |M|, ky = 0 we obtain bulk 
wavefunctions near the bottom of the conduction band: 


i’buik,tie-) cx 


/ cos fca;X-I-sin fc^^x \ /I 
i cos k^x — sin k^x 

o" 

0 


(A2) 


i’buik,iie) oc 


where kr 





0 


0 1 

cos kxX -1- sin k^x 

- 

1 

\—i cos kxX + |r^ sin kxX ^ 


\-*/ 

< A and /c « 

A 

|B| 

« A+ 


(A3) 


first terms in (A2|-(A3| correspond to the superposition 


Appendix B: Optical transitions in the case of 
strong electron-hole asymmetry 


In this section we analytically derive matrix elements 
of the optical transitions in the case of strong electron- 
hole asymmetry, i.e. \B — Vl B. For the simplicity 
we consider transitions between the Dirac point ky = 
0 and the bottom of the conduction band e — |AI| ^ 
\M\\/B'^ — jB. We also assume that typically ^ 
BM. For the edge states we can use the expression ([^, 

where A_ « A+ « > A_. 


2AB ’ + 

The bulk eigen-state of Hamiltonian for energy e > 
\M.\ is a sum of a right-moving wave, a left-moving wave 
and a term localized in the vicinity of the boundary: 


^{bulk) ^ ^ ^ ^^^d ( 3 ^) 

+ (B2) 

(b3) 

V''' = e-''"(^|Ei,s) + |iIi,s)), (B4) 


where k « A+, ^ k, A_. Amplitudes Ir 

and td can be obtained using the zero-boundary condition 
^(buik)^^ = 0) = 0. It yields 


tR = 


BkxK + iAi 
Bkj-K — iA4 


- 1 , 


td = - 


2BkxK 
Bka-K — iA4 


-2 


BikxK 

M 


(B5) 

























9 


Now the matrix elements can be calculated straight¬ 
forwardly. First, we find the overlap integrals defined in 
the section El 

_ ‘^ikx 1 _ 2Akx 1 , ^ 

A2 ^21;’ ^ ^ 

The matrix elements of the transitions can be calculated 
as 


ktP = |G| 


ktp = 


2\13-V\ 


A^k'^ 


\B-V\ 


A\B\X- 2\M\'^XlL^ \B\ 

1 2kl 


2A_ X^_Lr 


(B7) 

(B8) 


Appendix C: Matrix elements for the light incident 
at an arbitrary angle 

In this section we derive the matrix elements ([^ for 
light incident in the direction n^’^. We can take an aux¬ 
iliary orthonormal basis 

e®’'^ = cos (j) cos 9 + By sin 0 cos 9 — Bz sin 9 (Cl) 

(C2) 
i9 

(C3) 


= —Bj; sin (l) + By cos (p 

e®’'^ = n®’’^ = Bj; cos 4> sin d -|- sin p sin 9 + Bz cos 


The vector potential reads 

= e®’‘^Ao cos Wot — Sy’^Ao sinwot (C4) 

= Aq (cos (j) cos 9 cos ojA -\- sin p sin wot) (C5) 

Ay’‘^ = Aq (sin (j) cos 9 cos wot — cos p sin ujot) (C6) 

= “Alosin^coswot (C7) 

where the amplitude Aq = (cf. Eq. 4|. Since only 

in-plane components of the vector potentim yield the op¬ 
tical transitions, the case of arbitrary direction of light is 
equivalent to the case of elliptical polarization with sin 9 
standing for the eccentricity. 

The light-matter interaction Hamiltonian ([^ can be 
now rewritten as 


^e-A - 


Hie *‘^(cos0-|-l)-|- 
i7je*^(cos6l- 1) 


-f H.c. (C8) 


Time-reversal symmetry allows us to relate the matrix 
elements of Hi and hI 


‘Pedge 


H] 




,lk,; 


— \'4^edge, — s Hi 'Ipbulk,— 


(C9) 


This imply the relation between the matrix elements for 
different directions of the light 


e,4> 

= 


w, 


,e-*'^(cos 9 + 1) + n;l^e*'>(l - cos 9) 


(CIO) 


Appendix D: Derivation of kinetic equation 

In order to derive kinetic equation it is convenient 
to write a second-quantized Hamiltonian by expand¬ 
ing the field operator in eigen basis of Hbhz'- 4’ = 

E aky4ky + E [ck^+y4kl + ’ ^here is 

ky kx,ky 

the 4x1 wavefunction of the edge mode, 

4x1 wavefunctions in conduction and valence bands cor¬ 
respondingly. In this basis the Hamiltonian reads 


A 


Hbhz = ^s+ya's.k.As+y 

S,ky 


S,kx ,ky 


(c) At 


^^X,ky ^S,ky H“ 


E ' 

S,kx,ky 


(v) -t 

e) 1 V , , k k . 

kx.k-u S,kx,k^ 


He.—A — ^ ^ 


s,kx,ky 

W-siAz: + iAy) 

A) 


(Dl) 

(D2) 


K ky^e,kx+V 


where i. , el’ j, are the energies in the conduction 
and valence bands correspondingly. Here we disregard 
transitions between the edge states and the valence band 
as the frequency of illumination ojq is assumed to be much 
smaller than the difference between the Fermi level and 
the top of the valence band. 

The correlations of quantum fluctuations SA 
are given by the Green’s functions £>^^(1,1') = 

-i9{t -f)(^[sAM).sApir)]), D^^iiA') = i9{t’- 

t)([<5i„(l),<5i^(l')]), = -z({<5i„(l),M^(T)}). 

At zero temperature these Green’s functions are of the 
form 


D 


R(A) 


47r (nrUj)'^Sai3 - kakjd 
^ (n,w)2 ew2-|fc|2±i0 ’ 

= {D^ - £)"^)sign w 


(D3) 

(D4) 


However, due to the form of the Hamiltonian (D2) it is 
convenient to introduce axillary scalar fields 

Ag = Wg{Ax — iAy) + W—s{Az + iAy). 

Then the Green’s function for these fields i4f'(l,l') = 

-*d(t - 0)([^l(l),<5il(l')]), 4^(1,!') = lOit'- 
t)([5Agii),6AlA) 

take up the form 


= - 


*({ 5 X( 1 ),< 5 E(T)}) 


jjR{A,K) _ Dxx\'U}s + -|- Dyy\Ws — W—s|^~ 

i{Dxy +Dyx){wsW*_g-W-sW*g) (D5) 


We treat the interaction betw een electromagnetic field 
A and electrons given by (D2) as a perturbation, and 


use the Keldysh perturbation theory in order to derive 
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the retarded, advanced and Keldysh Green’s functions 
of electrons in the edge states. The Dyson equa¬ 
tions reads 


S{t -1') 


(D6) 


iidt - £s,ky)GXA^''^')~ 


f dt" + =0 

(D7) 

{-idf - es,ky) GXjt,t')- 

[ dt" + = 0 

(D8) 


The self-energy S can be represented as the sum of 
a classical contribution T,ind a contribution Egp due to 
quantum fluctuations <5.4, which is responsible for spon¬ 
taneous transitions 


E — ^ind + Egp, (D9) 

= (DIO) 


^R{A,K) 




2/')= ^(DfG^, + DfG«,) (Dll) 
V') = \ (Df G5 + DfGl, + 

(D12) 


Since the classical value of electromagnetic field (A) 
depends on time the self-energy Ei„d depend both on 
the sum of the times G = t + t' and their difference 
ta = t — t'. However, the dependence on the sum of the 
times describes the motion of electrons in high-frequency 
electromagnetic field, and we disregard it leaving only 
the dependence on the time difference, which is responsi¬ 
ble for induced transitions. After performing the Fourier 
transform over the difference time we obtain 


yR(A,K) , h \ — 
^ind.s — 


dTrW 

{nrUJof 


\Ws\^G^^J^'^\£ + UJ0, ky) + 


|rc_.pG5^-^)(e-a;o,fcp) (D13) 

The kinetic equation for the distribution function Usie) 
of electrons in the edge state may be derived in the stan¬ 
dard way by subtracting the Dyson equations with self¬ 
energy operators acting from the left (D7) and from the 
right (D8), integrating over ky, and using the relation 
G^ = (U^ - GA [1 - 2ns(e)]: 

2dtn{£) = i[l- 2n(e)] (E^(fcp) - ^^(ky)) - {ky) 

(D14) 


Using Eqs. (Dll) -(D13) and the expression for the 


Green’s functions in the conduction band G, 


R{A) 


c.s — 


1/ (e - £k^,ky ± *0), G^ = (G« - G^) [l - 2A,], where 
Ns (e) is the distribution function of electrons in the con¬ 
duction band, the summation over k^ in the self-energy 
operators is straightforward: 


^^nd{£^ky)-Y,fy^a{£,ky) = -2T:i '-'2 


(n^wo)^ 


kc,ky (e + Wo) 

(D15) 


^ind{£^ ^v) ~ 

IttW , 

- ^\ws\ i>s^ky{£ + i^o) [1 - 2A^(e-ka;o)] 

[UrUJo)^ 


(D16) 


Sfp(Gfcy) - ^7p{£,ky) = 

8i 


j (|wsp + |w_sp) i>c,ky{£ + w) N{e + uj)ujnr dtu 

(D17) 


Ssp(Gfcy) = 


8i 

“3 


J (|wsp + |w_sp) i>c,ky{£ + w) N{s + uj)ujnr duj 

(D18) 


The final kinetic equation can be obtained by substi¬ 
tuting the self-energy operators into (D14): 


dnJA 

dt 



nsje) - Nsje + oJo) ^^, ^ 

'£ind^s{£: £ T Wq) 

Ns{£ + u}) [1 - ns(e)] 

- 7 -\- 

Tsp{£ + U},e) 

„ 9~ / N IWsP 

VC^ky{£ + odo)-X2, 

^ nfujQ 

g 

- (llCsP -k |w_sp) U)l'c,ky{£ + 


(D19) 


(D20) 

(D21) 


Appendix E: EMF calculation 


From the kinetic equation (27), one can express the dis¬ 


tribution function fs and its angle-averaged value (/s)„ 


via the tunneling source Hg defined in (28): 


{fs)y, = 


ivkji + T 


-1 


^-1 


-n. 


lelvE 


dfo 

de 

(El) 


where k = (k||,A;z) is a parameter of spatial Fourier 
transform. 

The quasi-neutrality condition implies 


[ I \ 

J y *vk|| + j 


d( = Yl 

S 


n. 

ivk|[ -k 



d^ 

(E2) 
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Using the relation 


zvk|| +r 


and 


vE 


zvk|| +r ■ 


ik|| E 

'~kT~ 


1 - 


ip 'll Y 

nal equation for the electric field 


we obtain the fi- 


1kllE=^5: 


n. 


ivkii + T~ 




1 + kfl - 1 

(E3) 


The electric held averaged over the lateral direction x 
corresponds to the = 0 component: 


2D _ 

with cosw close to s- - — give the main contribution. 

After summation over spin index and integrating over 
^ = vf{\p\ — Pf) we obtain 

(Ey)^ (y) = 

7 r|Tp( 5 /r j" vkye^^^y 2sin ” 

L^PfVti tan ifo J v'^k^ cos'^ po + ^1 + k'^l'^ - 1 ^tt 

(E5) 


The EMF can be obtained as the integral of the mean 

Ly 

electric held £ = f {Ey)^ dy: 

-Ly/2 


{Ey}^ (y) = 


i 

2Lx 


E 


n. 


JkyV 


ivky cos (fi + T 
2sin^^l + fc2;2 

^1 + kill - 1 2^ 


(E4) 


where Us = 2f\T\'^5{x)5t (sp^.s “ [Is - «s] We 

perform averaging over y), assuming that the angles 


^ _ 2Fv\r\^5p r sin^ 2sin^^l + fcg/2 

tan ipo J v'^kl cos^ pQ + t-2 ^ 2^2 _i ^ 2 f 

(E6) 

The integration can be easily performed if Ly,A ^ It, 
i.e. when it is dominated by smal values of ky 

ri^- 1 ^ 1 ^ , A s, 
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